Removal of the Schwinger nonanaliticity in pair production 
by adiabatic switching of the electric field 
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Exact analytical solutions are presented for the time evolution of the density of pairs produced in 
the QED vacuum by a uniform electric field that is adiabatically switched on starting at — oo. Pair 
production is described by the Dirac-Heisenberg-Wigner function introduced before [Phys. Rev. D 
44, 1825 (1991)]. The explicit solution is obtained by an extension of the method of the spinorial 
decomposition to deal with a time- varying electric field. The main result of this work is that the 
pair density as a function of the field strength can be expanded into a power series. The essential 
singularity present in the Schwinger formula must be, therefore, attributed to the infinitely long 
duration of the process of pair creation by a time-independent field or alternatively by an abrupt 
switching of the field. 



I. INTRODUCTION 

The process of pair creation by the time independent 
and uniform electric field as described in our previous 
publication [Ij suffers from a serious drawback. The in- 
stantaneous switching on of the field leads to unphysi- 
cal transient effects clearly seen in Fig. 1 of In the 
present paper we solve analytically the problem of pair 
production by a uniform electric field that is switched on 
adiabatically. 

To explain what this means in physical terms, let us 
consider a large (boundary conditions are disregarded) 
capacitor charged so that the electric field is E. We shall 
compare two cases. In the first case, the electric field is 
primordial, it existed always. In the second case, the ca- 
pacitor has been continuously charged until aX t = Q the 
field reached the desired value E. The creation of pairs in 
the first case, i.e. for the time-independent electric field, 
was calculated by Schwinger [2!] whose analytical formula 
has a characteristic essential singularity as a function of 
the field strength S. 

The removal of the Schwinger singularity appearing in 
the title has the following meaning. Consider the ob- 
server who at time t — Q can measure the strength of the 
electric field and the number of pairs but does not have 
information about the past history. This observer can 
plot the number of pairs as a function of the field strength 
and find out whether this dependence shows the singular 
behavior characteristic for the Schwinger solution. We 
show with the help of exact analytic solutions that the 
singularity will not be found if the field has been turned 
on adiabatically. 

The formulas that describe the total pair production 
for an adiabatically switched field are analytic in the field 
strength — they can be expanded in the field strength into 
a convergent power series — for every finite value of the 
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parameter 6, that sets the time scale of the charging pro- 
cess. The Schwinger singularity in the field strength ap- 
pears now as a singularity in b. Thus, the essential singu- 
larity found by Schwinger is to be attributed entirely to 
rather unphysical assumptions that the field has infinite 
duration or that it has been suddenly switched on. 

The production of pairs by an adiabatic perturbation 
may seem impossible since there is an energy gap of 2mc^ 
between the vacuum state and the state with one pair. 
Therefore, it would seem that under an adiabatic per- 
turbation the vacuum will not change. However, these 
doubts should disappear when one realizes that pair pro- 
duction by a constant electric field is due to tunneling. 
Arbitrarily weak electric field can produce an arbitrarily 
large difference in potential energy for a sufficiently large 
separation. When this difference in potential energy ex- 
ceeds the value 2mc^, pairs are created. 

The problem of pair production will be formulated in 
terms of a field-theoretic generalization of the standard, 
nonrelativistic, Wigner function. The generalization was 
introduced in [j] and called the Dirac-Heisenberg-Wigner 
(DHW) function. This formulation was further devel- 
oped in [i-[l2. It allows for a description of pair creation 
as a time evolution of the vacuum state. The advantage 
of using the DHW function is its applicability to the ex- 
pectation values of field operators and not only to the 
one-particle wave functions that are solutions of the Dirac 
equation. However, the analytic solutions of the one- 
particle Dirac equation are still very useful. The method 
of spinorial decomposition introduced in enables one 
to construct from these solutions the components of the 
DHW function. We shall describe two exactly soluble 
examples. 

In the first example the field is switched on at a con- 
stant rate b — the time dependence is governed by a simple 
exponential. In the second example the field is switched 
on at a decreasing rate — the time dependence is gov- 
erned by the hyperbolic secant squared. In both cases 
the signature of the emerging Schwinger singularity is 
the shrinking of the radius of convergence to zero when 
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b tends to zero. We also show that in the first case the 
power expansion of the exact solution coincides with the 
straightforward perturbation theory. 

The evolution equation for the DHW function has 
many solutions, most of them do not represent physi- 
cal states. The DHW function shares this property with 
the standard Wigner function in nonrelativistic quantum 
mechanics — not every function over the phase-space is a 
Wigner function. For example, the Wigner function can- 
not be concentrated too much in position and momen- 
tum since this would contradict the uncertainty relations. 
There is, however, a class of problems for which one can 
be sure that the Wigner function represents a physically 
admissible state. These are the problems of time evolu- 
tion with the initial value being a genuine Wigner func- 
tion. The DHW function that is the solution of the evo- 
lution equation will correspond to a physical state — the 
one that evolved from the initial state according to the 
Schrodinger equation. This is precisely the method that 
is used in this work. To describe the time evolution un- 
der the influence of adiabatically switched on fields we 
shall always start with the easily computable exact DHW 
function — the one that describes the free vacuum state. 

In nonrelativistic quantum mechanics the Wigner func- 
tion is defined as a Fourier integral of a product of wave 
functions and in many cases this integral can be explicitly 
evaluated. In contrast, the calculation of the DHW func- 
tion in quantum field theory requires the summation of 
contributions from an infinite number of states of virtual 
particles in the vacuum. The only case where this can 
be easily done is the vacuum state of the free field. For 
that reason in our analysis of the time evolution we shall 
choose the DHW function of the vacuum as the initial 
value. 

In general, the evolution equations for the components 
of the DHW function, even in the approximation where 
all radiative corrections are neglected, form a set of 16 
integro-differential equations and their analytic solution 
is a hopeless task. However, in the case of a spatially 
uniform electromagnetic field the evolution equations be- 
come tractable since they reduce to linear partial differ- 
ential equations. These equations allow for the applica- 
tion of the method of characteristics leading to ordinary 
differential equations. In turn, these ordinary differential 
equations in certain cases can be solved analytically by 
spinorial decomposition introduced in Ref. [l|. 

In the present paper we construct a complete set of so- 
lutions of the evolution equations for the DHW function 
describing the evolution of the vacuum in the presence of 
a uniform but time-dependent electric field. With the use 
of these solutions we find the evolution when the initial 
state is the free-field vacuum. Our solution clarifies the 
mechanism responsible for pair production and may ex- 
plain some details of the vacuum decay studied recently. 

In Sec. HI] we recall basic definitions and we derive 
a simplified form of the evolution equations best suited 
for our analysis. In Sec. IIIII we exhibit the underly- 
ing 0(3) symmetry that enables us to express the solu- 



tions of the evolution equations in terms of spinors. The 
method of spinorial decomposition is extended in this 
section to cover a time-dependent field. In Sec. lIVl using 
this method, we find an analytic solution of the spinorial 
equations when the electric field is turned on exponen- 
tially. In Sec. lVIII we construct the analytic solution from 
the solutions of the Dirac equation found by Sauter 13;]. 
Finally, in Sec. I Villi we compare the time evolution of 
the QED vacuum in both cases underscoring the main 
features. 



II. EVOLUTION OF THE DHW FUNCTION 
IN A TIME- VARYING ELECTRIC FIELD 

The DHW function is defined as follows B: 
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*„(r + r,/2,t),#t (r-r,/2,t) $), (1) 



where the square bracket denotes the commutator of the 
quantized Dirac field operator ^ with its hermitian con- 
jugate and the phase ip is the line integral of the vector 
potential 
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- / dXr]-A{r + Xr],t). 
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The physical meaning of this phase factor is best seen 
by rewriting Eq. ([T]) with the use of the gauge covariant 
shift operator 



ai3{r,p,t) = J (fir] exp{-ip-q/h) 

e^/'-^^^{r,t),e-^/'-^i'l {r,t) 



X ($ 



(3) 



where D = 'V — ie/hA{r, t). The equivalence of this form 
of W and the previous one follows from the relations 



ei^/^-D = cxp 



i- / dXri-A(r±Xri,t) 



(4) 



The phase factor has been introduced in 3] to secure 
gauge invariance of the DHW function. That is why only 
the electromagnetic fields and not the potentials appear 
in the evolution equations. This also means that the 
argument p of the DHW function is the kinetic not the 
canonical momentum. 

The most general DHW function has 16 real compo- 
nents. The 4x4 matrix Wa/s can be decomposed into 16 
matrices Q. We found it convenient to use the matrices 
introduced originally by Dirac [T3 |. 

3 3 

W{r,p,t) = j[fo + ^Pif,. + (T-ga + ^Pi(T-gi]. (5) 
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In the special case when the initial state is the vacuum 
state and the applied electric field is spatially uniform, 
only the following quantities do not vanish: the mass 
density f3{r,p,t), the spin density go{r,p,t), the elec- 
tric current density gi{r,p, t), and the magnetic moment 
density g2ir,p,t). Moreover, as shown in P, 0], only 
three independent coefhcients {wi,W2,W3) are needed to 
express all these quantities, namely 



n 



(6) 

where n is the unit vector in the direction of the electric 
field, p± is the component of p perpendicular to the field, 
and E± is the transverse energy, E± = ^/m'^c^ + c^p^. 
The factor -1/2 has been introduced to have the vector 
(?i'i,W2,W3) normalized to 1. 

The evolution equations for the three coefficients were 
derived in [ij and they are: 
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These equations differ from Eqs. (11) of Ref. [l| only by 
the time-dependence of the electric field. We introduced 
here the dimensionless variables, namely 

ehcE{t) 



Ej_t 



El 



(8) 



where p\\ is the component of p in the direction of the 
electric field. 

We shall expose the group-theoretic content of the evo- 
lution equations for the vector w{q,T) by writing them 
with the use of the spin-one matrices S in the Cartesian 
basis, 



{dr + £{T)d^) w{q, r) = 2iS-n{q)w{q, r), (9) 



where 



niq) = {l,0,q), S-n{q) 



g 

-q 1 
0-10 



(10) 



In this way we may interpret the evolution as a rotation 
in the three-dimensional vector space p^ . 

In order to convert partial differential equations into 
ordinary differential equations we use an extension of the 
method of characteristics to systems of equations. How- 
ever, this time the conversion goes in opposite direction 
as compared to our previous work. Namely, we shall elim- 
inate the derivative with respect to q, not the time deriva- 
tive. This is achieved by replacing the vector w(q,T) by 
a new vector X(q, r) with the momentum variable q dis- 
placed by the vector potential A{t), 



q^q + A{T), A{t) 



dr'Sir'). (11) 



Formally, this is accomplished by the following substitu- 
tion: 



w{q, r) = exp [A{T)dq] X{q, r) 



(12) 



When this expression is inserted into Eqs. ([9]) we find 
that the components of this new vector satisfy ordinary 
differential equations with q playing the role of a fixed 
parameter. 



2iS-n{q-A{T))X{q,T). (13) 

This equation was obtained with the use of the relation 

exp [-A{T)dq] niq) exp [Air)dg] = ?i(q - A{t)). (14) 

The solution of the initial value problem will proceed 
by construction of the fundamental matrix J^(q, r) for 
the set of equations ([13]). This matrix is built (cf., for 
example, |1Q] ) from a set of three linearly independent 
solutions (Afi(q,r), Af2((?,r), A'3(g,T)) of Eqs. 
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(15) 



The fundamental matrix normalized to the unit matrix 
at To is 



J^N{q,T,TQ) = J^{q,T)T ^{q,To). 



(16) 



The solution of the initial value problem for the vector X 
with an arbitrary initial condition X{q, tq) set at r = tq 
is therefore 



= TNiq,T,To)X{q,To). 



(17) 



Our aim is to find the time evolution of the vacuum 
state. The DHW function for this state is built from the 
following set of Wi functions [3] : 
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W2ip)^0, W3{p) 
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(18) 



where Ep — ^ rn?-c^ -\- c^p^ is the total energy. Thus, in 
our dimensionless units the vector describing the vacuum 
has the form: 



y(g) = -^(i,o,<z), 

Eq 



(19) 



where Eq — \J\ + q^ represents the dimensionless total 
energy. Choosing Af(g,To) = V{q) as the initial state we 
obtain the time evolved DHW function by inserting the 
formula dTT]) into Eq. (HI]) 



w{q, t) =pN{q + A{t),t, To)V[q + A{t)). 



(20) 



All we need now to solve the initial value problem is the 
fundamental matrix. In the next section we show how 
to construct this matrix using the method of spinorial 
decomposition. 
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III. SPINORIAL DECOMPOSITION 

The set of three differential equations leads to third 
order equations for each component studied in Ref. [6]. 
As we have shown in [l|, one may replace third order 
equations with second order equations by employing the 
method of spinorial decomposition that exploits the 0(3) 
group symmetry. We shall extend this method here to 
cover the case of the time- varying field. To this end, 
we introduce the following representation for the compo- 
nents of the vector X in terms of a two-component spinor 
^ and the Pauli matrices, 



X"^ = *tcr^*. (21) 



One may check by inspection that the evolution equations 
for the vector X follow from the evolution equations 
for the spinor ^ and its conjugate. 



^*(9,t) 



t) = -z^Hq, T)cT-n{q - A{t)). 
These two coupled equations 



-^^J-{q,T) = i{q- A{T))^/j_{q,T) +i^j^{q,T), 
^'0+(9,t) -i{q- A{T))%l;+{q,T) +iip_{q,T) 



(22a) 
(22b) 

(23a) 
(23b) 



can be transformed to the following two separate second 
order equations for the upper ?/;_ and lower compo- 
nents of the spinor. 



+ l + iq~AiT)fT^£ir) 



^^iq,T)=0. (24) 



These equations coincide with the second order equations 
obtained from the Dirac equation in the time-dependent 
potential A{t) along the z axis after the separation of 
the spatial variables through the following substitution 



^/'(r, t) — exp(ip • r)ip{t). 



(25) 



Therefore, known analytic solutions of the reduced Dirac 
equation ([M)) will lead to analytic solutions for the 
Wigner function. Linearly independent solutions of 
Eqs. (IM)) serve as building blocks in the construction 
of the spinors from which we shall construct the funda- 
mental matrix (|16p . The most convenient choice of two 
spinorial solutions is when they are orthonormal. Let us 
suppose that we found one solution ^'i (q, r) of Eqs. (j22al) . 



*i(9,t) = 



/ 

-g 



(26) 



The minus sign in the lower component is chosen for fu- 
ture convenience. Then, the second solution of Eqs. (|22ap 
can always be chosen in the form 



«'2(g,T) 



(27) 



It is automatically orthogonal to the first solution and it 
is normalized if the first solution was normalized. 

The following formulas, introduced in show how 
to produce three orthonormal vectors Xk from two or- 
thonormal spinors and 4*2: 

Xi = -i (*Icr*2 + "flfT^i) , (28a) 
X2 = -i (-i^\(T^2 + i*J(T$i) , (28b) 



(28c) 



The orthonormality of these vectors is a result of the 
following identity valid for every set of four spinors. 



= 2 (4-1*4) (^3^2) - (^1*2) (*J«'4 



(29) 



The fundamental matrix built from orthonormal vectors 
is orthogonal. Therefore, in the formula ([T5|) we 
may replace the inverse matrix by the transposed matrix. 
We have chosen the minus sign in the definition I^E^ to 
have the vectors {X i, X2, X^) forming a right-handed 
set, Xi - {X2 X Xs) = 1. Therefore, the determinant of 
the fundamental matrix is equal to 1 {T describes a pure 
rotation). 

We may now use the three orthonormal vectors 
to build the matrix T. The components of this matrix 
[TtI are real and imaginary parts of various quadratic 
expressions of the two functions / and g, 



-m\]g\ 23[/5] gg* - ff* 



(30) 



We need a specific form of the function A{t) to proceed 
any further. This will be done in next sections. 



IV. ELECTRIC FIELD SWITCHED ON AT A 
CONSTANT RATE 

The adiabatic switching on of the electric field will be 
described by an exponential prefactor. Thus, the electric 
field will be given by the formula S{t) = £ exp(6T)n with 
6 > 0. It follows from the Maxwell equations that this 
field is produced by the current b£{T). For sufficiently 
small values of b this current can be neglected in com- 
parison with the current associated with the produced 
pairs. 

Assuming that the initial time tq lies in the remote 
past (tq = —00), we obtain 



A{t) = -£:6-iexp(6r). 



(31) 



This choice of the time dependence results in the follow- 
ing second order equations (j24p for the upper and lower 
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spinor components: 

[^3 + 1 + (9 + r T *f e"^] {q, r) = 0. (32) 

These equations are solvable in terms of confluent hyper- 
geometric functions. The two independent solutions for 
■(/)_ have the form 



exp(zTi?q)e 



'is/2 



V2jE^„ - qE, 



(33a) 



:iFi(l-i^^-^,l-2i^;is] , (33b) 



/(9,t) = 



s/2^El + qEg 

.E, 
b 



where s = 2£b~'^e^'^. The lower components are, in gen- 
eral, linear combinations of both f*{q,T) and g{q,T). 
However, for our choice p3p of the upper components 
the lower components "0+ contain only one term and the 
two orthonormal solutions of Eqs. (|22p are (p6)) and (l27t . 
This is the simplest pair of solutions because each com- 
ponent contains only one confluent hypergeometric func- 
tion. Orthogonality of the two spinors is self-evident but 
the check of normalization requires the use of the evolu- 
tion equations ()22p . Since these equations preserve scalar 
products of spinors, it is sufficient to check the normal- 
ization at any value of r and we choose t = —oo or s = 0. 
At this value, the confluent hypergeometric functions re- 
duce to 1 and the normalization |/p + |gp = 1 follows 
immediately. It may be worth noticing that the normal- 
ization condition for an arbitrary time generates some 
quadratic identities for the confluent functions that seem 
to be unknown. The matrix J-^ will be used in Sec. IVlto 
determine from Eq. (|20p the time evolution of the vacuum 
state. 

To elucidate the connection between the results ob- 
tained for the time-independent field in [l[ and the time- 
dependent field in the present case, we concoct one "su- 
per" equation that encompasses these two cases. The 
evolution equation (j22ap will have the form: 



dT 



*(g,r) = icr-n [q + £ 



,exp(fer) - exp(6ro) 



*(9,t), 
(34) 



where tq specifies the initial moment at which 5* describes 
the field-free vacuum. The exact solution of this equation 
can still be found. All we have to do is to replace q by 
q — exp(6ro)/6 in ([33]). There are two mutually exclusive 
limiting procedures that lead either to the solution of 
Eqs. (22) in [l| or to the time-dependent solution given 
by (|33p. The first limit is when tq = and 6 = and 
it leads to a solution singular at £ = 0. This solution 



reproduces the Schwinger singularity as shown in ^ . The 
second limit is when tq — >■ — oo. This solution is analytic 
in £ but it has a singularity when 6 — > 0. Thus, there is 
no direct connection between these two limiting cases. 



V. TIME EVOLUTION OF THE QED VACUUM 

One must be careful with taking the limit tq — > — oo in 
Eq. ((20)) because of the oscillatory terms present in the 
fundamental matrix even in the limit of vanishing field. 
However, these terms cancel out when the initial state is 
the vacuum, 



lim T'{q,To)V{q) 

T0-!--OO 



(0,0,-1). 



(35) 



The final expression for the vector Wv{q,T) that evolved 
from the vacuum at tq = — oo is obtained by combining 
the formulas pO | . (|30)) . and (|35|) . 



w^{q,T) = T{q + A{t),t){0, 0,-1) 

= (-23?(/r),25(/r),i/i 



' \~9\' 



(36) 



where / = /(g -I- yl(r), r) and g = g{q + A{t),t). This 
vector satisfies the evolution equations ([7]) and the vac- 
uum initial condition. The validity of these statements 
follows from our general procedure but it can also be 
checked by a direct (albeit tedious) calculation. 

The density of pairs n(g, r) produced by the electric 
field [ll] was shown in Q to be 



n{q,T) = 1 - 0{q,T), 



(37) 



where 0{q, r) is the overlap between the DHW vector at 
time T given by (I36p and the initial vacuum vector (|19p. 



Oiq,T)^Viq)-w^{q,T) 

9(i/T - m - mfr 



E„ 



(38) 



Since both vectors V{q) and w^{q,T) are normalized, 
1 — 0{q, t) is simply one-half of the squared difference 
between the initial and the evolved vectors. The rela- 
tion (137p has been guessed from numerical calculations in 
Ref. [3| . In the Appendix we give a derivation of this im- 
portant formula and show that it is valid for an arbitrary 
electromagnetic field. In the present case /o vanishes be- 
cause the charge density in a state that evolved from the 
vacuum under the action of a uniform field vanishes at 
all times. Therefore, the general formula (|A8p reduces to 



mf3{r,p,t) +p-gi{r,p,t) 



2Er, 



(39) 



With the use of Eqs. ^ and (jlSp the second term can 
be identified with (minus) the overlap function and the 
formula (I37p is thus proven. 
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FIG. 1. The pair density n{q, r) for the exponentially growing 
field at time t — and for b — 0.05 plotted as a function of 
the dimensionless momentum q at the values 0.1, 0.2, 0.3, and 
0.4 of the dimensionless field strength £. 



The number of pairs should not change under the si- 
multaneous reversal of the field and the parallel momen- 
tum because it is only the relative orientation of these two 
vectors that matters. This property directly follows from 
the fact that under the change of signs £ — >■ —£, q ^ —q 
the functions / and g are interchanged / — ?■ —g, g — > — /. 



In Fig. [T] we show the pair density as a function of q for 
different values of the field strength £. In contrast to the 
results of Ref. there are no transients. The number 
of pairs is a smooth function of q. 

The most striking result is the absence of an essential 
singularity in the dependence of the number of produced 
pairs on the field strength £. All functions appearing in 
the formula as long as 6 7^ 0, can be expanded into 
a convergent power series in £ because the confluent hy- 
pergeometric function 1-^1(0, c; z) is an analytic function 
of its three arguments. However, it has poles in c at all 
negative integers —n (including 0). Since in the formulas 
psp the parameter c is a function of £, 

c=l±2.V^i±IlEW, (40) 

the radius of convergence R of the series in £ can be 
determined from the conditions c — and it reads 

R = by/l+lfTWJl. (41) 

Thus, the perturbative expansion is convergent only for 
£ smaller than the R. In the limiting case, when 6 — )■ 0, 
the radius of convergence goes to zero as it must because 
of the essential singularity that reappears in this limit. 

The expansion in £ for t = gives the following for- 
mula up to the sixth order in £: 



O(,,0) = l ^£^ ^3 



E^{A + b^) £:9(4-f 62)2(1 + ^2) £;12(4 + 62)3(1+ 62)2(4 ^_ 952) 

_f 5 . . — . _f ^ (42) 



E\^(\ + 62)4(1 + ^2)3(4 + 962)2(1 + 4^,2) £:18(4 + 52)5(1 + ^2)4(4 + 9^2)3(1 + 462)2(4 + 2562) 

where 6 — h/ Eq and 

A^2 = 1/2, (43a) 

N^ = q{fb + P), (43b) 

Ni = 3/8^ - 448q2 + 112 + P{nmq^ - 104) + 6''(2164g2 _ 313) + h^{4,12q^ - 106) + 6**(36(?^ - 9)) , (43c) 

iV5 = (7/2(6(-126720(72 + 65472) - b^{ll0b28q^ - 180048) + 6^(755440(?2 - 39900) + 6^(1332348?^ - 425793) 

-I- 69(7951369^ - 379041) + 6"(193144g2 - 124383) + 6"(24804q2 _ 17631) + h^^{12mq^ - 972)), (43d) 

iVe = 1/16(48660480(7'' - 50282496g2 + 3559424 + 62(-1305722880(?'' + 924008448(?^ - 30885888) 

+ 6'*(-11614668800g'* -f 11574500352g2 _ 593933568) -I- 66(-24204609536g'* + 41551564800g2 - 2746675968) 
+ 6^(21767894400^'' + 52683929664?^ - 5392492368) + 6i°(151238728256q'* - 21806265408q^ - 3779458248) 
+ 612(236416765448(7'* - 131339980524^^ + 2556452997) + 6i'*(183464805216g'* - 150072654336^^ + 6447749532) 
+ 61^(78987949264(7"* - 83787726072g2 + 4744001262) + 6*^(198257727049'' - 25548039744?^ + 1733216492) 
-f 6^" (3204676008?^ - 4441116924?^ + 339299253) -I- 6^2(293855040?" - 426785760?^ + 34399080) 
+ 624(11664000?'* - 17496000?2 + 1458000)) (43e) 
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FIG. 2. The analogous plot as in Fig. 1 obtained from the 
sixth-order perturbative expansion. 



In Fig. 2 we show that the six terms of this series repro- 
duce reasonably well the exact results. 

The structure of the denominators in Eq. ([^^ is a di- 
rect consequence of the series expansion of the confluent 
hypergeometric function 



iFi(a, c;z) = 1 + - — + 
c 1! 



a{a + 1) 



c(c+ 1) 2! 



a{a + l)(a + 2) 
' c{c+l)(c + 2) 31 



+ 



(44) 



In our case c is given by Eq. PO)) and the expansion of 
the expressions l/{n + c) into powers of £ gives 



b{{n + l)bT2i) 2iq{{n+l)bT'2if 



± 



1)262)2 



(45) 



Substituting here the consecutive values of n we repro- 
duce the factors in the denominators of (l42l). 



Even though each term of the expansion into powers 
of £ is finite, for b = the series becomes only an asymp- 
totic one. In order to make this point absolutely clear, we 
calculated for q = the first 8 terms of the perturbation 
series for the number of produced pairs for two values of 
b. 



_ 1 2 63 4 25513 g 43884005549 g 235474684928440509 

~^ " ~ 10 ~ 2600 9802000 ^ 7834150480000 ^ 26610807325949600000 

1776480040832261689003353 12 189228367844501729409089588397558027 14 

^ 180781584001131645584000000 ~ 19653421122106203692350489578400000000 
171936287478376823113691752600797604830985578228669 . ^g. 

19571206822106176531328917776928027793134400000000000 + ' ^ ■> 

^ ^_^c2 '^^ cA 869 g 334477 8 59697183 34429291905 12 14631594576045 ^4 

"^128 ^ 1024 ^ 32768 ^ 262144 4194304 ^ 33554432 

68787420596367165 ^ig 
+ 2147483648 ^ ) 



These expansions leave no doubt that the first series is 
convergent for £ <1 (the coefficients tend to zero) while 
the second is only asymptotic (the coefficients grow ex- 
ponentially). The transition from the convergent series 



to the asymptotic series is the mechanism for the reap- 
pearance of the Schwinger singularity when the adiabatic 
switching is removed. 
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VI. PERTURBATION THEORY 



The solution of the evolution equation p3)) for the vec- 
tor X{q, t) in perturbation theory is most easily obtained 
by expanding this vector in the eigenbasis of the ma- 
trix appearing in ([7]) ■ This matrix has eigenvalues and 
±2iEq. Its orthonormal eigenvectors are: 



e±(g) 



1 



(48a) 



(48b) 



We shall seek the solutions of equations ([T^ in the form 

X{q,T) = xo{q,T)eQ{q) 

+ x+{q,T)e+{q)+x^{q,T)e^{q). (49) 

The evolution equations for the expansion coefficients are 





dr 








dr 




-(9,r) 



2iEqX+{q,T) 

iA{t) (2qx+{q,T) - V2a;o(g,r)) , (50b) 

-2i£'gX+(q,r) 
iA{T)(2qx-{q,T)-V2xo{q,T)j , (50c) 

where A{t) — A{T)/Eq. In order to prepare these equa- 
tions for the perturbative solution we shall multiply the 
second equation by exp{—2iEqT) and the third equa- 
tion by exp(2z£'qT) and then integrate all three equation 
with respect to r from — oo to r. Assuming that ini- 
tially at T = — cxD the system is in the vacuum state, i.e. 
xo{q, ^oo) = 1, x±{q, — oo) — 0, we obtain after integra- 
tion by parts the following set of three coupled integral 
equations: 

xo{q,T) = 1 - iV2 dri 

J — oo 

X A{ti) {x^{q,Ti) - x+{q,Ti)) , (51a) 
x+{q,T) = -I rdTie^'^'^"-"^) 

J — OO 

X A{tx) (2qa;+(q,ri) - %/2 a;o(g, n)) , (51b) 

J —OO 

X A{ti) (2qx^{q,Ti) ~ V2xo{q,Ti)'^ . (51c) 

The solution of these equations by iteration produces the 
perturbative series in powers of the field strength which 



coincides with the power expansion of the analytic solu- 
tion. We checked that perturbation theory reproduces 
the formula (H^ . Note that the function A{t) blows up 
for b — 0. Hence, the integral equations (15(1)) become 
invalid which is another manifestation of the Schwinger 
singularity in the absence of adiabatic switching. 



VII. ELECTRIC FIELD SWITCHED ON AT A 
DECREASING RATE 

To complete our study we shall show now that the 
main conclusion of our investigation — the removal of the 
essential singularity — holds also when the time depen- 
dence of the electric field is different. To this end we 
use analytic solutions of the Dirac equation in a uniform 
time-dependent field discovered by Sauter [l^ and ana- 
lyzed in more detail by Sommerfeld [l9| . The time depen- 
dence of the electric field in the Sauter solution is £{t) — 
£sech^{bT/2)n so that A{t) = -25(1 + tanh(6r/2)). 
In this case, as in the previous one, the electric field 
is turned on gradually, starting from the zero value at 
t = — oo. 

The second order equations ([24|l for the Sauter field 
have the following form: 

r d^ ^ , 2f (1 +tanh(6r/2)),p 



T^ . M ^{q,T)=0. (52) 
cosh (or) 



Following Sauter and Sommerfeld we replace the time 
variable by u = (1 -I- tanh(5r/2))/2 and write down the 
solutions in terms of hypergeometric functions. Two or- 
thonormal solutions (|26l) and (p7)) of the first order equa- 
tions ([25)1 are constructed this time from the following 
two solutions of the second order equation for -0- : 



M*"(l - u) 



^/2JeI ~ qEq 



'X2F1 { 4i— + ia - - 4i— + ia - 1/3,1 + 2ia; u ) , 



(53a) 



g*{q,u) = - 



u-*"(l - u) 



03 



V2jE^q + qEq 



f £ £ \ 

X2F1 I 4i— — ia + — 4i— — ia -|- 1 — 2ia;u j , 

(53b) 

where 



The number of pairs constructed according to Eqs. (j37p 
and (1381) is an analytic function of £ but this time the 
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FIG. 3. The pair density n{q) for the Sauter field at time 
T — and for 6 = 0.05 plotted as a function of the dimension- 
less momentum q at the values 0.1, 0.2, 0.3, and 0.4 of the 
dimensionless field strength £. 




-20 -15 -10 -5 

T 

FIG. 4. The field strength £ shown as a function of time. The 
lower curve corresponds to the exponential dependence and 
the upper curve to the square of the hyperbolic secans. 

radius of convergence is four times smaller than the value 
(|4T|) for the exponential field. For the Sauter field we re- 
produce the formula for the number of produced pairs at 
t = -|-oo that were obtained by Narozhnyi and Nikishov 

^ sinh {^-n{a- ^)) sinh + ^ (a ~ ^)) ^ ^^^^ 
sinh {2na) sinh (27r/3) ' 

where /? differs from /3 by having the sign of £ reversed. 



VIII. COMMON PROPERTIES OF THE TWO 
SOLUTIONS 

The two exactly solvable cases differ considerably but 
they share common features confirming the main message 
of our work. We have chosen the parameters in such a 
way that at r = the electric field in both cases attains 




-20 -15 -10 -5 

T 

FIG. 5. The pair density n{q, r) for the exponential field 
(lower curve) and the Sauter field (upper curve) plotted as a 
function of time. Since at r = both fields reach the same 
value, the number of pairs is almost the same. 

the same value £ and the rate of growth of the field in the 
remote past is b. In this way we will be able to compare 
directly the results for the two choices of £{t). In Fig. |3] 
we show the time dependence of the electric field in the 
two cases. 

There is a direct correspondence between the functions 
that solve the Dirac equation. The hypergeometric func- 
tions ((33)) are related to the confluent hypergeometric 
functions by the following limiting procedure 

T — > — oo, £ — > (X), %£e^'^ — const = h^s. (56) 

Since at r = the electric field attains in both cases 
the same value f , we expect that the number of produced 
pairs will not differ significantly. Indeed, the results of 
the numerical calculations shown in Figs. [5] and [5] fully 
confirm this expectation. The number of produced pairs 
is slightly higher for the Sauter field because, as shown in 
Fig. m this field dominates over the exponentially grow- 
ing field. 

The analytic character of the number of pairs function 
is seen in Fig. [6] where we plotted the dependence on £ 
in the neighborhood of zero. The parabolic character of 
both curves clearly shows that the expansion in powers of 
£ begins in both cases with the quadratic term. The dif- 
ference between the two cases diminishes with decreasing 
b. 



IX. CONCLUSIONS 

We have shown with the help of exact analytic solu- 
tions that adiabatic switching of the electric field removes 
the essential singularity which is present in the Schwinger 
formula for pair creation. This observation opened the 
way to a perturbative treatment which gives the results 
coinciding with the power expansion of the exact solu- 
tion. The method of spinorial decomposition that en- 



10 




0.04 



0.02 



FIG. 6. The pair density n{£) for the exponentially growing 
field at the time t — for b = 0.1 and q = plotted as a 
function of the dimensionless field strength £ for the expo- 
nential field (lower parabola) and for the Sauter field (upper 
parabola) . For comparison we also show a typical nonanalytic 
behavior representing the Schwinger singularity. 



abled us to find exact solutions is quite general and might 
be applicable to a class of similar problems. 
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Appendix A 

The starting point is the expansion of the Dirac field 
operator (cf., for example, ^22i] ) into creation and anni- 
hilation operators (c=l), 




-ip-r /h 



(Al) 



We allowed for an arbitrary time dependence of the cre- 
ation and annihilation operators, so that the validity of 
this formula is not restricted to free fields. Next, we in- 
vert the Fourier transform and use the properties of the 
spinors u and v to obtain 

, / m I f (Pr 

b{p,s,t) = J—uT{p,s)- 



(27r;i)3/2 



(A2) 



With the use of the completeness relation for the spinors 
Ufip, sX(p, s) = — '-, (A3) 



expressed in the original Dirac notation, we obtain the 
following formula for the number of electrons with mo- 
mentum p: 



N{p,t)^Y.^\p,s,t)b{p,s,t) 



i'{r + r]/2,t). (A4) 



P3m + pi(T-p + Ep 



2E„ 



Now, we replace the product of field operators by the 
commutator and anticommutator and use the canonical 
equal-time anticommutation relations 



1 r 
2 



T,/2),*J,(r- 



rj/2) . (A5) 



We have not assumed that the field operators obey the 
free Dirac equation since we allowed in Eq. (IA1[) for an 
arbitrary dependence of the creation and annihilation op- 
erators on t. However, in the presence of an external 
field we must introduce modifications required by gauge 
invariance. Following the procedure that led to the for- 
mula ([3]) we replace the shift operators in Eq. (jA4|) by 
their gauge invariant counterparts @ and extract the 
phase factor e~^^ . The gauge invariant number operator 
Ng obtained in this way is 



Ng{p,t)=J2b\p.s,mp,s,t) 



■^{r + r]/2,t). (A6) 



{2TTh) 

Pzm + pi(T-p + Ep 



2E„ 



The expectation value of Ng{p,t) can then be expressed 
in terms of the DHW function. The quantity of interest 
is the pair density in phase space n(r,p,t) which in the 
present case is equal to the density of electrons. There- 
fore, in the formula (jA6p we may drop the integration 
over r and the normalization of the phase-space volume 
element {2Trh)^ to obtain the density of pairs per unit 
phase-space volume in the form 



n{r,p, t) 

Tr{(/93m- 



pi(T-p + Ep)W{r,p,t)} 



2 En 



(A7) 



Due to the orthogonality of the Dirac matrices under the 
trace only the components /3,gi, and fo may contribute 
and we obtain the final formula 



n{r,p,t) 

mf3{r,p,t) +p-gi{r,p,t) 



= 1 



Epfoir,p,t) 



2Eri 



(A8) 
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The number of antiparticles calculated in the same man- 
ner is 

n{r,p,t) 

^ ^ _^ mf3ir,p,t)+p-giir,p,t) - Epfo{r,p,t) ^^^^ 

2Ep 

The difference between the two expressions is the charge 
density 

nir,P,t) ~ n{r,p,t) = fo{r,p,t), (AlO) 
as was to be expected. 
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